﻿Modelare şi simulare Seminar 3 SEMINARUL 3 3 Procese Markov În cadrul proceselor aleatoare un loc deosebit îl ocupă procesele Markov caracterizate prin aceea că apariţia unei anumite stări este condiţionată doar de un număr determinat de stări anterioare Dacă numărul acestor stări anterioare este r, atunci este vorba de un proces Markov de ordinul r Procesele Markov de ordinul 1 ocupă un loc important în studiul proceselor de trafic ce caracterizează reţeaua de telecomunicaţii în ansamblul ei Înseamnă că dacă un sistem se află la … şi dacă el se comportă ca un momentele de timp ,, , ,tt t……în stări 01j01(), (), , (),jNt Nt Nt… sistem Markov de ordinul 1, atunci probabilitatea ca la un moment următor t starea lui să fie 1j+ (Nt) este determinată doar de starea lui din momentul anterior, adică: 1j+ == = =…= 10011{( ) () , () , , () }jjPNt kNt n Nt n Nt n+ (3 1) ==== 1{( ) () } ()jjkPNt kNt n p t+ Cu alte cuvinte aceasta se exprimă prin formularea: starea viitoare k depinde doar de starea prezentă n sau starea prezentă înglobează în ea întregul trecut Se mai obişnuieşte să se spună că procesele Markov sunt procese fără memorie Pentru orice sistem este caracteristic vectorul de stare ( tp El este structurat ca un vector linie, ) având ca elemente probabilităţile absolute ( rilor ale sistemului (n0n≤), ≤∞ ) npt ale tuturor stă adică: (3 2) p…… () 01((),(),, (),)ntptpt pt= cu 0( ) 1pt≤≤ şi n() 1npt= ()n∑ Dacă se consideră evoluţia unui astfel de sistem între două momente successive t şi tt+, Δ înseamnă că o stare ktΔ oarecare poate fi obţinută printr-un salt executat în intervalul dinspre oricare altă stare a sistemului fie ea 0,1,2, , sau n 1Modelare şi simulare Seminar 3 În figura 3 1 se prezintă grafic această evoluţie corespunzătoare unui proces Markov de ordinul 1 geţi sunt Stările sistemului numerotate începând cu 0, sunt reprezentate prin cercuri, iar prin să figurate eventualele tranziţii între stări Pe fiecare săgeată sunt notate probabilităţile de tranziţie: ,kkp ,nkp nk()tt+Δ ()t 0,kp 2,kp 1,kp 0 2 1 ()t()t ()t Fig 3 1 Tranziţii între stări la modelul Markov =+Δ= = (3 3) ,() {( ) () }nkptPNttkNtnΔ Cu alte cuvinte, pentru ca la momentul tt+Δ sistemul să se afle într-o stare nouă k, trebuie ca plecând din starea prezentă nnk→ la momentul t să aibe loc tranziţia: , adică: starea ( ) {starea ( ) salt }kt t nt n k+Δ=∪ ∩ → ()n(3 4) ceea ce în termini probabilistici se poate scrie ca: = ⋅ Δ ∑ ,() ()(kn)nkptt ptp t+Δ (3 5) ()n Această relaţie reprezintă modul de exprimare matematică a caracterului Markovian al procesului, în care starea viitoare kn nu depinde decât de starea prezentă Relaţia 3 5 poate fi scrisă şi sub formă matriceală Toate probabilităţile de tranziţie p sunt ,nk grupate într-o matrice de tranziţ ii T (vezi tabelul 3 1), o matrice pătratică ce are ca dimensiune 2Modelare şi simulare Seminar 3 numărul total de stări posibile ale sistemului În această matrice suma elementelor fiecărei linii este egală cu unitatea Se obţine astfel relaţia matriceală: t ( ) () () (0)tt t sau t+Δ = ⋅= ⋅ppTppT (3 6) care face posibilă determinarea vectorului de stare a sistemului în orice moment de timp, pe baza rii vide matricei de tranziţii şi a vectorului de stă 3 1 Procesul general de naştere şi moarte Orice reţea de telecomunicaţii funcţionează conform unor procese în echilibru statistic, în care apelurile (cererile de serviciu) se nasc, aşteaptă servirea, trăiesc o viaţă atingând o anumită vârstă, pentru ca apoi în final să moară Procesele de naştere şi moarte ca subclasă importantă a proceselor Markov de ordinul I se adaptează perfect modului de apaiţie şi de servire a traficului în sistemul global de telecomunicaţii Aceste procese sunt caracterizate de un spaţiu discret de stări şi de un parametru temporal continuu Prin definiţie, un proces de naştere şi moarte autorizează pentru un sistem, care la momentul t se '( ' )tt t> n=1kn=+ia la momentul următor numai către stările adiacente află în starea k, tranziţ k't sau 1kn=− şi evident că starea poate fi menţinută şi la momentul următor Evoluţia dinamică a unui sistem conform unui proces de naştere şi moarte poate fi reprezentată ca în de stări fig 3 2 sub forma unui lanţului Markov liniar cu număr finit 1N+ p p p 1,NN− 0,1 1,nnp− ,1nn+ N 1n−1N− n1n+ 0 1 1,0p p ,1p− NN ,1nnp−1,nn+ Fig 3 2 Lanţ Markov liniar (proces general de naştere şi moarte) ia de tip nn→ Pentru simplificarea desenului nu s-a reprezentat decât pentru starea n tranziţ (fluxul intern de probabilităţi p); de asemenea nu este marcat nici intervalul 'tttΔ=în care se ,nk − a pot produce tranziţiile, controlate prin termeni de tipul ()pttΔ ,nk De fapt ,()nkpt apare ce o rată măsoară chiar probabilitatea tranziţiei nk→ schimbării de stare şi care multiplicată prin tΔ 3Modelare şi simulare Seminar 3 Pe durata intervalului tΔ nu posate avea loc decât cel mult o singură tranziţie, iar aceasta este permisă doar între stări adiacente Deci conform relaţiei 3 5, se poate scrie că: (3 7) +Δ =Δ +Δ + − Δ − Δ 11, 11,,1 ,1( ) () () () () () n nnnnnnnnn nnpttptptptptptptpt−− +++ − În regim staţionar, probabilităţile stărilor sunt independente de timp, iar valoarea probabilităţii de tranziţie depinde, în mod direct proporţional, de lungimea intervalului în care are loc tranziţia şi de rata p a schimbă reală pozitivă, adică: ,nkrii de stare, care este o constantă Δ+ Δ (3 8) ,,() ()nknkptptotΔ= Cu ()otΔeste notat indicele Landau definit astfel: ()otΔ (3 9) lim 0= 0tΔ→ tΔ Aşa după cum s-a recizat anterior, pentru orice proces Markov se poate forma matricea T de tranziţie, ale cărei elemente sunt probabilităţile de tranziţie între stări În condiţiile echilibrului ia echilibrului global scrisă ca în relaţia (3 10) şi care probabilistic exprimă statistic se obţine ecuaţ faptul că suma fluxurilor de ieşire de probabilităţi dintr-o stare este identic egală cu suma fluxurilor de intrare: += + (3 10) ,1 ,1 1 1, 1 1,()nnn nn n nn n npp p p p p p+−−−+n+ N1 1n− n1n+ 1N− n-1 0 1∑ 4∑ 3∑ 2∑ Fig 3 3 Suprafeţele de conservare a fluxurilor de probabilităţi (I) 4Modelare şi simulare Seminar 3 p 1,nnp− ,1nn+ ∑ 1n+ n-1 n p 1,nn+ ,1nnp− Fig 3 4 Suprafeţele de conservare a fluxurilor de probabilităţi (II) iile echilibrelor locale după cum urmează: Se pot deduce deasemenea ecuaţ == 11,0 0 0,11 1,, 1;;NNN NNNpp pp p p p p−−− (3 11) == 11, ,1 ,1 11,;;nnn nnn nnn nnnpp pp pp pp−−− + ++ Ultimele două relaţii din 3 11 indică faptul că în ecuaţia echilibrului global termenii din cei doi membrii sunt egali doi câte doi, adică pentru o stare tranziţ niile ascendentă şi descendentă sunt echiprobabile Dacă se consideră variaţia continuă a parametrului timp, atunci plecând de la (3 5) este adevărat că: − +Δ − ,,()( ) 1()()nkkkkkpt p tpt t ptΔ limlim ( )( ) limpt pt=⋅+ ∑ 00 0nktt tΔ→ Δ→Δ→ tΔΔΔ (3 12) ()nktt≠ Se formează în consecinţă matricea infinitezimală Q, ale cărei elemente sunt constante reale şi anume: ,()nkptΔ ,limnkq= 0tΔ→ tΔ (3 13) − ,()1kkptΔ ,limkkq= 0tΔ→ tΔ Relaţia (3 12) devine: ' +q⋅ ∑ ,,() () ()knnkkkkptptqpt=⋅ (3 14) ()nk≠ ia viitorului: Sau sub forma cunoscută ca şi ecuaţ ' Q() ()tt=pp (3 15) 5Modelare şi simulare Seminar 3 Trebuie precizat că elementele matricei infinitezimale Q (vezi tabelul 3 1) se deduc uş or din cele ale matricei de tranziţii T 1 QTI (3 16) ()=− tΔ 3 2 Procesul de naştere pură Un proces de naştere pură are o evoluţie ca cea prezentă în figura 3 5 în care s-a considerat un proces trunchiat, adică s-a presupus un set finit de 1N+ stări −λ Δ 1t− λΔ 1t− 1t−λ−λΔ 1 Δ 11nt− 11nt+ nλΔ1N− 01t−λ 1 Δ 0 n-1 1n− N1 n1n+ 1N− tλΔ Δ tλ 1t− nΔ nλ 1N− 0tλΔ t Fig 3 5 Lanţul tranziţiilor într-un proces de naştere pură În raport cu o stare ( )Nt n= , pentru un moment de timp ttΔ, în intervalul imediat următor sunt posibile doar următoarele două evoluţii: ¾ Apare o naştere: +Δ = + = = Δ (3 17) ,1{( ) 1 () } ( )nnprob N t t n N t n p t+ ¾ Nu se întâmplă nimic: = = = Δ (3 18) ,{( ) () } ( )nnprob N t t n N t n p t+Δ Ceea ce înseamnă că: Δ+ Δ= (3 19) ,1,() ()1nnnnptpt+ Procesul este omogen şi probabilităţile de tranziţie nu sunt funcţii de timp Ele sunt proporţionale cu tΔ, adică: 6Modelare şi simulare Seminar 3 Δ= Δ+ Δλ ,1() ()nnnpttot+ (3 20) tλ − Δ+ Δ ,()1 ()nnnpttoΔ= reprezintă rata medie de naşteri adică numărul mediu de cereri de serviciu ce apar în unitatea de nλ timp Mărimea intervalului de tranziţie tΔ este în aşa fel aleasă încât probabilitatea de a avea mai mult de tΔ o naştere în tΔ să fie un infinit mic de ordin superior lui , adică ()otΔ Tabelul 3 1 Matricea T pentru procesul de naştere pură ri StăStări în tt+Δ în t 0 1 n-1 n n+1 N-1 N 0 1t−Δ tΔλ 0 0 0 0 0 0λ0 1 0 1t−Δ 0 0 0 0 0 1λ n-1 0 0 1λ −tΔλ 0 0 0 1n−1n− n 0 0 0 1tλ ΔtλΔ 0 0 n−n n+1 0 0 0 0 1tλ 0 0 −Δ 1n+ N-1 0 0 0 0 0 1tλ −ΔtλΔ 1N−1N− N 0 0 0 0 0 0 1 Tabelul 3 2 Matricea Q pentru procesul de naştere pură ri StăStări în tt+Δ în t 0 1 n-1 n n+1 N-1 N 0 − λ 0 0 0 0 0 0λ0 1 0 − 0 0 0 0 0 1λ n-1 0 0 λ − λ 0 0 0 1n−1n− n 0 0 0 λ λ 0 0 n−n 7Modelare şi simulare Seminar 3 n+1 0 0 0 0 λ − 0 0 1n+ N-1 0 0 0 0 0 λ − λ 1N−1N− N 0 0 0 0 0 0 1 in relaţiile: Plecând de la ecuţiile viitorului şi folosind elementele din Q se obţ ' λ 000() ()pt pt=− ' =−λλ t (3 21) 11() () ()nnn nnptptp−− ' =λ 11() ()NNNpt p t−− Procesul de sosire adecvat cererilor de serviciu se caracterizează prin independenţa faţă de starea sistemului a ratei de naştere λ că: n, ceea ce înseamnă λ pentru orice nN∈ (3 22) ,= nλ Presupunem că la momentul iniţial 0t =, starea sistemului este = Înseamnă că: 0(0) 0Nn= ⎧ 01, dacann= (3 23) (0)p= n⎨ 00, dacann≠ ⎩ ie Poisson: Se poate deduce uşor că distribuţia căutată este o distribuţ n ()tλ t −λ ()pte= n (3 24) !n Ea precizează de altfel probabilitatea condiţionată ca n naş teri să apară în intervalul (,)0, t dacă starea iniţială este caracterizată de o populaţie nulă, adică 0n= 0 Aplicaţia 1 Linia unui abonat telefonic se poate găsi în două stări: liber sau ocupat, notate cu 0, respectiv 1 Ştiind că duratele celor două tipuri de stări urmează distribuţii de probabilitate exponenţiale de medie 154α =, respectiv 15β=, să se determine probabilitatea de stare ()pt şi ()pt în 01 i (0) 0, 7p= funcţie de probabilităţile iniţiale de stare: (0) 0,3p= 0 ş1 8Modelare şi simulare Seminar 3 Rezolvare: , ⎧ =− ()()() 11nnn nnptptptλλ −− ⎪ ,⎨ =− () () () ⎪ 101pt pt ptαβ ⎩ , =− + t () () () 11nnnnnptptpμμ ++ , =− + () 001()()pt p t ptαβ ⇒ sist de ecuaţii diferenţiale a viitorului: ' ⎧ =− + () ⎪ 00() ()1ptptptαβ '⎨ =− ⎪ 101() () ()ptptptαβ ⎩ - relaţia de normare: 01() () 1pt pt+= ' =− + − t () 00() (1 ())0ptpt pαβ ' β ++ = () 00()()pt ptαβ β ++ = 00() ( ) ()sP s P sαβ ββ 0()Ps== SSαβ ++++ αβ Cu soluţia: β ()tαβ−+ () 0pt Ce=+ + αβ 145 1αβ +=+== 55 5 1 1β 5 0, 2== = + 15αβ β =−= 0, 3 0, 3 0, 2 0,1C=− + αβ Cu soluţia: Penru t =0 avem condiţia: (0) 0 3p= 0 ββ 0, 30, 3CC=+⇒=− ++ αβαβ 9Modelare şi simulare Seminar 3 t− + 0() 0,2 0,1pte⇒= ' =− t 10() () ()1ptptpαβ 01() 1 ()pt pt=− ' =− −pt () () 11() 1()1ptptαβ ' α++ = () 11() ()pt ptαβ α tαβ−+ () 1()pt Ke=+ + βα pentru t = 0 condiţia iniţială 1(0) 0, 7p= αα +⇒=− 11(0)(0)pKKp⇒= βαβα++ 4 4α 5 = 0,8== 145βα + + 55 =− 0, 7 0,8 0,1K=− t− () 10,8* 0,1pt e= Definiţia 1 Traficul este raportul dintre durata medie de serviciu şi durata medie a interapelului, ceea ce înseamnă că e o cantitate fără dimensiuni: 1secμ () A= 1seλ c () Unitatea care tinde să se generaliz este ERLANG-ul Definiţia 2 Traficul oferit sau traficul de intrare este nr mediu de cereri ce se prezintă în sistem, cu rata de , în intervalul de apelare egal cu timpul medie de serviciu 1μ apelare λ cereri secunda 1 secundeAλμ=⋅ ()() Dar cererile nu pot exista în afara surselor cele produc Aceeaşi definiţie poate căpăta forma: Traficul oferit este nr mediu de surse simultan ocupate în decursul unui interval cât durata medie de , dacă rata lor de apelare este λ servire 1μ surse ocupate secunda 1 secundeAλμ=⋅ ()() 10Modelare şi simulare Seminar 3 Definiţia 3 Orice cerere acceptată de sistem are la dispoziţie un circuit sau dispozitiv de prelucrare, pe care îl ocupă pe toată durata lui de prezenţă în sistem, ceea ce înseamnă că: Traficul scurs (prelucrat) este nr de resurse simultan ocupate pe durata timpului mediu de serviciu 1μ pentru prelucrarea cererilor oferite cu rata λ resurse ocupate secunda 1 secundeAλμ=⋅ ()() Aplicaţia 2 În figura de mai jos este prezentată diagrama de ocupare a 3 resurse, funcţionând într-un grup comun observate în decurs de 20 ore Se cere: a grupului a ) Diagrama de ocupare totală de angajare b ) Durata totală c ) Intensitatea traficului scurs b ) Durata totală de angajare: 3 σσ=++ ∑ 123iσσ 1i= =++++= 13,520,52,519,5oreσ =+ + + = 23 3,5 2,5 1,5 10,5oreσ =+ + +++= 321,52,532112oreσ 11Modelare şi simulare Seminar 3 c ) Intensitatea traficului scurs: R 11As iiσσ σσ==⋅+⋅+⋅= ()() ( ) ( ) ∑ 1iTT T= 1 [1 1 2 0, 5 3 0, 5 2 1 3 1 2 1 1 2 2 3 1, 5 2 0, 5= ⋅+⋅ +⋅ +⋅+⋅+⋅++⋅+⋅ +⋅ + 20 32ore 1 2 0, 5 0, 5 1 2 1 3 0, 5 2 0, 5 0, 5 0, 5 2] 1, 6E+⋅ + ++⋅+⋅ +⋅ + + + = = 20 = () ∑∑ iiiσσ Nr sursei 1 2 3 Valori medii Durata totală 9,5 10,5 12 Intensitate trafic:1,6 de angajare σ i Nr de 5 4 3 angajărin i Rata de 5/20 = 0,25 4/20 = 0,2 3/20 = 0,15 Rata medie de angajare apelare:0,6 individuală λ i Durata totală 9,5/5 = 1,9 10,5/4 = 2,26 12/3 = 4 Durata medie de de servicii tsserviciu:2,66h/apel i 'iniσ =; ts= iλ i Tn i Aplicaţia 3 În urma observării a 4 resurse pe durata a 6 minute s-au trasat diagramele din figura de mai jos, în care: 1 - ocupat; 0 – starea de neutilizare a fiecărei resurse Se cere: a ) Diagrama de ocupare totală a grupului b ) Durata totală de angajare c ) Intensitatea traficului scurs 12Modelare şi simulare Seminar 3 Stabiliţi intensitatea medie a traficului servit de respectivul grup pe parcursul intervalului considerat Rezolvare: 1− i n σ (durata în div ˆA ⎡⎤ iτ ii miniλ i ⎣⎦ de servicii) 1 2 2,5 2/6 = 1/3 2,5/2 = 5/4 2,5/6 = 5/12 2 3 3 3/6 =1/2 3/3 = 1 3/6 =1/2 3 2 3,5 2/6 = 1/3 3,5/2 = 7/4 3,5/6 = 7/12 4 4 4,5 4/6 = 2/3 4,5/4 = 9/8 4,5/6 = 9/12 =3/4 Total 11 13,5 11/6 27/22 – durata 9/4 medie de serviciu 44 1 1 113 5 27 9A s ===+++== ( ) ( ) (2 5 3 3 5 4 5)ii iσσ = ∑∑ 124 1166iiTt T== 9 E 2 25As== ∑ i= 4 i iσ () iiσ () 1 0,5 0,5 2 2 4 3 3 9 4 0 0 Total 13,5 minute Aplicaţia 4 1 Accesul unui grup de terminale la reţeaua de telecomunicaţii se face prin intermediul unui mod de comunicaţie care dispune de o singură linie de date de o capacitate de 10 Mbiţi/sec Dacă traficul oferit de terminale spre a fi transferat în reţea este de: ( i ) 600 pachete / sec cu câte 1000 de biţi pachetul ( ii) 8 Mbiţi/sec 13Modelare şi simulare Seminar 3 exprimările în unităţi Erlang: /*(mesaje oră Lbiţiλ) () B= i ( / sec) *3600Dbiţ În cazul nostru: 600 / sec*1000 6pachete = () 50, 06iBE== 10*10 100 6 i 8*10 / secbiţ () 60,8ii BE== i 10*10 / secbiţ Notă: Trafic scurs după definiţia 3 RR ==i () ∑∑ siAiσσ 11ii== 1 Intensitatea traficului după definiţia 4 () ∑ Kiiσ 11siA= == ∑ iσ 1iTTT= 2 Timpul mediu de servicii caracteristic iσ *ts= i in 3 Nr de angajări în unitatea de timp (rate de angajare) ,in = iλ T 14